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computed  as  applied  to  the  “representative”  phosphorus  smoke  for  growth  toward  equilibrium 
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smoke  is  established.  Numerical  values  indicate  that  yield  factors  are  highly  dependent  on  the 
environmental  relative  humidity. 
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PREDICTING  THE  DROPLET  SIZE  AND  YIELD  FACTORS  OF  A PHOSPHORUS 
SMOKE  AS  A FUNCTION  OF  DROPLET  COMPOSITION  AND  AMBIENT 
RELATIVE  HUMIDITY  UNDER  TACTICAL  CONDITIONS 


I.  INTRODUCTION. 

The  development  of  a particle  size  prediction  model  satisfies  two  basic  informational 
needs  of  aerosol  physicists:  first,  the  prediction  model  provides  a means  for  characterizing  the 
physical  properties  of  particles  comprising  an  aerosol;  secondly,  the  model  illustrates  possible 
avenues  toward  the  aim  of  artificial  particle  size  control. 

The  goal  of  this  study  is  threefold:  (l)To  develop  a predictive  relationship  between 
the  size  of  a phosphoric  acid  droplet  and  the  ambient  relative  humidity.  (2)  To  explore  the 
kinetics  of  condensational  growth  toward  stable  equilibrium  and  to  investigate  the  validity  of  the 
assumption  that  vapor  diffusional  equilibrium  exists  between  the  environment  and  the  surface  of 
the  phosphoric  acid  droplet  under  tactical  conditions.  (3)  To  formulate  a predictive  relationship 
between  the  yield  factors  of  a phosphorus  smoke  and  the  relative  humidity. 

In  order  to  realistically  model  the  growth  and  predict  the  droplet  size  of  an  aerosol 
such  as  phosphorus,  several  physiochemical  processes  must  be  investigated.  However,  due  to  the 
hygroscopic  nature  of  the  phosphorus  oxides  and  their  tendency  to  generate  solution  droplets, 
one  such  process,  condensation-evaporation,  assumes  a dominant  role  in  the  determination  of  the 
droplet  size  under  tactical  conditions. 

The  analysis  of  solution  droplets  has  received  a considerable  amount  of  attention 
from  such  researchers  as  Byers,1  Low,2  and  Friedlander.3  To  effectively  analyze  the  relationship 
between  vapor  condensation  and  the  size  of  a phosphoric  acid  droplet,  a combination  of 
equilibrium  thermodynamics  and  aerosol  kinetics  is  employed.  With  data  for  the  water  activity, 
surface  tension,  and  density  of  phosphoric  acid  solutions,  the  droplet  (phosphoric  acid) 
saturation  ratio  is  formulated  from  basic  thermodynamic  considerations.  By  applying  the  Kelvin 
and  solution  effect  to  a phosphoric  acid  droplet,  the  droplet  saturation  ratio  is  explicitly 
expressed  as  a function  of  droplet  diameter  and  acid  content.  To  incorporate  the  solution  effect, 
an  empirical  relationship  between  the  water  activity  and  water  mole  fraction  of  a phosphoric 
acid  solution  is  formulated.  Vapor  diffusional  equilibrium  is  assumed  between  the  surface  of  the 
phosphoric  acid  droplet  and  the  environment  whereupon  a predictive  relationship  between 
relative  humidity  and  droplet  diameter  follows. 

To  verify  the  validity  of  vapor  diffusional  equilibrium,  the  kinetics  of  condensational 
growth  is  explored.  Employing  the  modified  growth  law  developed  by  Carstens,4  a quasi-steady 
state  analysis  based  upon  a balance  between  heat  conduction  and  vapor  condensation,  the  growth 
of  a phosphoric  acid  nucleus  can  be  delineated.  The  growth  law  is  modified  to  apply  to  the 
growth  regimes  dominated  by  molecular  and  continuum  phenomenon.  From  the  modified  growth 
law,  relaxation  times  (times  for  growth  toward  stable  equilibrium)  are  constructed.  Numerical 
results  indicate  that  phosphoric  acid  droplets  “track”  with  the  environment. 

Subsequently,  yield  factors  for  a phosphorus  smoke  are  displayed. 
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II.  RELATIONSHIP  BETWEEN  THE  SATURATION  RATIO  AND  DIAMETER  OF  A 

SOLUTION  DROPLET. 


In  order  to  establish  the  relationship  between  the  relative  humidity  and  droplet 
diameter,  the  vapor  (water)  saturation  ratio  over  a droplet  containing  a nonvolatile,  soluble, 
dissociating  solute  must  be  developed.  Theoretically,  if  the  saturation  ratio  of  the  environment  is 
equal  to  1,  that  is, 


where 


P = P, 


so 


(1) 


Se  - environmental  saturation  ratio,  alternatively  referred  to  as  the  ambient 
relative  humidity 

P = actual  vapor  pressure  of  the  environment 

Pso  = equilibrium  vapor  pressure  determined  by  the  temperature  of  the  environment, 
vapor  pressure  over  a plane  surface  of  pure  water 

then  liquid-phase  water  should  coexist  in  equilibrium  with  gas-phase  water.  However,  nucleation 
experiments  indicate  that  water  droplets  will  not  form  until  the  saturation  ratio  of  the 
environment  exceeds  3 (P  > 3 Pso),  an  extremely  supersaturated  state.* 


The  reason  for  this  discrepancy  between  experiment  and  theory  can  be  explained  in 
terms  of  the  Kelvin  effect.  The  Kelvin  relation  relates  the  vapor  pressure  over  a curved  surface  to 
the  vapor  pressure  over  a plane  surface  of  the  same  liquid  and  concludes  that  the  vapor  pressure 
over  a curved  surface  is  greater  than  the  vapor  pressure  over  a plane  surface.  From  a physical 
viewpoint,  the  vapor  pressure  over  a surface  of  liquid  is  determined  by  the  energy  necessary  to 
separate  a volatile  molecule  on  the  surface  from  the  attractive  forces  exerted  by  its  nearest 
neighbors  located  in  the  layer  adjacent  to  the  surface. ^ The  attractive  number  ratio,  defined  by  the 
ratio  of  the  number  of  molecules  in  the  layer  adjacent  to  the  surface  (Ns_j)  to  the  number  of 
molecules  in  the  surface  layer  (Ns),  is  smaller  for  the  curved  surface  as  compared  to  the  plane 
surface.  Consequently,  the  net  energy  of  attraction  is  effectively  reduced  over  the  curved  surface 
thereby  augmenting  the  escaping  tendency  of  the  volatile  molecules. 


Mathematically,  this  relationship  between  vapor  pressures  may  be  given  as:1,2 


Ln 


4a,M1 

RTp'd 


(2) 


*This  is  true  only  in  the  absence  of  foreign  material. 
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where 


Pj  = vapor  pressure  over  the  curved  surface  characterized  by  diameter,  d 
Ps  = vapor  pressure  over  a plane  surface 


a'  = surface  tension  of  the  surface  with  units  ergs/cm^ 
p ' = density  of  the  solution  droplet 
Mj  = molecular  weight  of  water 
R = the  universal  gas  constant 
T = temperature  of  the  system 

The  Kelvin  law  states  that  the  natural  logarithm  of  the  ratio  of  the  two  vapor  pressures  is 
inversely  proportional  to  the  diameter  of  the  droplet;  all  other  terms  are  considered  constant.*  By 
transforming  the  Kelvin  law  into  its  exponential  equivalent,  it  is  clearly  seen  that  P^  > Ps  for  all 
diameters,  d.  The  implications  are  that  when  the  environmental  vapor  pressure  is  just  able  to 
sustain  a plane  surface  of  liquid  (Se  = 1),  it  is  not  yet  capable  of  sustaining  a droplet  with  a 
curved  surface.  Only  when  the  environment  is  supersaturated  (P  > 3 Pso)  can  droplets  exist  in 
equilibrium.  If  this  were  the  only  phenomenon  controlling  the  vapor  pressure  over  a droplet 
containing  a soluble  nonvolatile  solute,  then  the  existence  of  a droplet  in  an  urisaturated 
environment  would  be  a thermodynamically  impossible  state.  However,  much  air-quality 
degradation  in  unsaturated  environments  is  due  to  droplet  formation  (haze)  exemplified  by 
sulfuric  acid  droplets. 

The  phenomenon  that  allows  droplets  to  exist  in  equilibrium  with  an  unsaturated 
environment  is  commonly  referred  to  as  the  solute  effect.  The  solute  effect  relates  the  vapor 
pressure  over  a plane  surface  of  liquid  with  soluble,  nonvolatile  solute  to  the  vapor  pressure  over 
a plane  surface  of  liquid  without  solute.  Because  of  the  solute  effect,  the  vapor  pressure  over  a 
liquid  with  solute  is  less  than  the  vapor  pressure  over  a liquid  without  solute.  The  physical 
reasoning  that  supports  this  conclusion  is  as  follows: 

Vapor  pressure  is  primarily  a colligative  property**  dependent  on  the 
relative  number  of  water  molecules  available  for  entry  into  the  gas 
phase.  With  the  addition  of  solute  molecules,  the  relative  number  of 
water  molecules  is  effectively  reduced;  therefore,  a reduction  in  vapor 
pressure  consequently  follows.3 


In  this  study,  variations  in  o and  p result  in  second-order  variations  in  saturation  ratios. 

** 

This  is  strictly  true  for  ideal  solutions  only. 
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In  ideal  solutions,  the  solute  effect  can  be  mathematically  modelled  by  Raoult’s  law, 


where 


Ps  *lPso 


(3) 


Ps  = vapor  pressure  over  a plane  surface  containing 
liquid  with  soluble,  nonvolatile  solute 

Xj  = mole  fraction  of  solvent  (1) 

Ps0  = vapor  pressure  over  the  pure  liquid  with  plane  surface 

With  the  addition  of  solute  molecules,  the  solvent  mole  fraction,  Xj,  is  reduced 
from  its  maximum  value  of  unity  to  some  fractional  value,  thereby  resulting  in  a proportionate 
reduction  in  the  vapor  pressure  of  the  solution  as  inferred  from  equation  (3).  From  past 
experimental  work,  data  on  the  vapor  pressure  over  a plane-surfaced  solution  of  phosphoric  acid 
and  water  has  been  amalgamated  for  a wide  range  of  acid  concentrations  at  25°C5’6  The  data 
indicate  a greater  vapor  pressure  reduction  due  to  the  addition  of  a given  amount  of  solute  as 
compared  to  the  value  predicted  by  Raoult’s  law  and  support  the  conclusion  that  phosphoric 
acid  solutions  are  nonideal.  The  nonideality  of  the  solutions  may  stem  from  two  basic  physical 
mechanisms:  (l)the  dissociation  of  H3PO4  in  aqueous  solution  and  a subsequent  increase  in  the 
relative  number  of  solute  molecules  and  (2)  molecular  dipole-dipole  interactions  between  H3PO4 
and  H^O. 


Nevertheless,  the  inappropriateness  of  Raoult’s  law  as  a model  for  phosphoric  acid 
solutions  is  quite  evident  and  obviates  its  application  in  this  study. 

A more  cogent  model  for  the  solute  effect  in  nonideal  solutions  can  be 
mathematically  represented  by 


Ps  ^left  Pso 

where  Xjeff,  the  effective  solvent  mole  fraction  or  water  activity,7  accounts  for  the  nonideality 
of  the  solution. 

Through  a joint  effort  between  Mr.  Theodore  Tarnove*  and  this  author,  the 
functional  relationship  between  X]eff  and  Xj  has  been  mathematically  formulated  for  phosphoric 
acid  solutions.  In  order  to  account  for  the  total  domain  (0<  Xj  < 1),  two  straight-line  fits  were 
applied,  resulting  in  the  following  relationship: 


Xleff  " 

2X,  - 1 

1 > 

X1 

> 

.577 

Xlet'f  = 

.22X, 

0< 

X1 

< 

.577 

*Tarnove,  Theordore,  Munitions  Division,  Chemical  Systems  Laboratory,  Aberdeen  Proving  Ground, 
Maryland  21010,  February  1978,  private  communication. 
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The  general  expression  for  the  saturation  ratio  over  the  droplet  containing  soluble 
dissociative  nonvolatile  solute  is  obtained  by  the  substitution  of  equation  (4)  into  equation  (2), 
resulting  in  the  following  thermodynamic  relationship: 


In 


Pd  4o'Mj 
P^~  RTp'd 


- In 


xleff 


(5) 


III.  SATURATION  RATIO  OVER  PHOSPHORIC  ACID  DROPLETS. 

For  phosphorus  smoke,  the  empirical  relationships  for  Xjeff  may  be  employed,  along 
with  the  substitution  X,  = 7- — , where  n,  and  n9  are  the  number  of  moles  of  solvent  and 

n2+ni  1 z 

solute,  respectively,  thereby  resulting  in  the  multivalued  functional  form  for  the  saturation  ratio 
over  the  surface  of  phosphoric  acid  droplet,  where: 


In 


d _ 


so 


4a'Mj 

RTp'd 


- In 


1 + 


2m 


nl  " n2 


1 > Xj  > .577 


(6) 


Pd  4a'Mj 

In  = In 

Pso  RTp'd 


7rd3 

With  the  additional  substitution, = n2  ^2  + nl  ^1  >*  which  states  that  the  total  volume  of  the 

6 

droplet  is  associated  with  the  volume  of  the  solute  and  the  volume  of  the  solvent,  then 
equations  (6)  and  (7)  become  with  some  mathematical  simplification: 


n2 

4.55  + 4.55  — 


*1  J 


.577  < X,  > 0 


(7) 


Pd  _ 4a'Mj 

In  = - In 

Pso  RTdp 


1 + 


1 2n2  V j 

?rd3  - 6n2(V2  + Vj) 


1 > Xj  > .577 


(8) 


In 


_Pd 

Pso 


4a 'Mi 

- In 

RTdp' 


4.55 


27.3n2  Vj 
7rd3  - 6n2  V2 


.577  > Xj  > 0 


(9) 


*From  density-data  analysis,  V2  (the  mole  volume  of  H3PO4)  and  Vj  (the  mole  volume  of  H20)  are  reasonably 
constant  over  all  degrees  of  acid-concentrated  solutions. 
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The  above  particular  expressions  relate  the  saturation  ratio  over  a phosphorie  aeid  droplet, 

pi 

i.e.,  Sa  = , to  the  fundamental  properties  of  that  droplet,  i.e.,  the  diameter,  d and  the  number  of 

PSO 

moles  of  aeid  n2  contained  in  the  droplet.  The  first  term  on  the  right  (positive)  represents  a vapor 
pressure  increase  due  to  the  Kelvin  effect;  and  the  seeond  term,  a subtraetive  expression,  represents 
the  vapor  pressure  reduction  due  to  the  solute  effeet. 

Trd^ 

With  the  approximation  that  the  droplets  contain  dilute  solutions,  i.e.,  — g-  » n2V2, 

then  the  logarithmic  term  may  be  expanded  in  a Taylor  series  expansion  retaining  the  first  term 
only:  ln(  1 +X)  = -X  where  X<<1.  This  circumstance  would  be  applicable  only  to  equation  (8) 
where  the  mole  fraction  of  solvent  approaches  unity  and,  therefore,  with  this  approximation 
equation  (8)  becomes: 


Pj  _ 4a'M}  1 2n  2 V t 

P^"  RTp'd 


(10) 


This  expression  is  referred  to  as  the  Kohler  expression  and  is  valid  only  for  dilute  solutions.3  While 
it  is  not  a good  model  for  highly  concentrated  aeid  droplets,  as  in  a phosphorus  smoke,  some 
valuable  information  may  be  gleaned  from  this  expression.  In  particular,  the  competition  between 
effect  (the  first  term  on  the  right)  and  the  solute  effect  is  elearly  revealed  as  a function  of  diameter. 
For  “small”  diameters,  the  solute  effeet  dominates  the  Kelvin  effect,  resulting  in  saturation  ratios 
less  than  1 (see  figure  1).  As  the  diameter  increases  and  the  solution  becomes  dilute,  the  two  effects 
become  equal  in  magnitude  (Sj  = 1)  and  thereafter  the  Kelvin  effect  dominates,  resulting  in  a 
droplet  saturation  ratio  greater  than  1.  As  the  diameter  of  the  drop  increases  to  infinity,  both 
effects  become  negligible,  physically  representing  a plane  surface  of  pure  liquid,  whose  saturation 
ratio  equals  1. 


To  effectively  model  phosphorus  smoke,  consisting  of  concentrated  phosphoric  acid 
droplets,  equations  (8)  and  (9)  will  be  employed  for  analysis. 

To  gain  additional  insight  from  equations  (8)  and  (9),  a graphical  representation  is 
displayed  on  a Sj  versus  diameter  axes  system  for  constant  mole  of  solute.  The  result  is  a family  of 
curves,  eaeh  eurve  constrained  by  a fixed  amount  of  solute.*  This  procedure  has  been  applied  to  a 
phosphorus  smoke  where  the  moles  of  aeid  contained  in  each  droplet  has  been  determined  in  the 
following  manner  (see  figure  1 and  table  1). 

Using  a cascade  impaetor,  measurement  of  aeid  droplet  sizes  at  known  relative 
humidity  is  performed;  and  with  knowledge  of  the  relationship  between  aeid  concentration  of 
the  solution  and  relative  humidity,  a determination  of  number  of  moles  of  solute  contained  in 


*In  general  a two-component  system  has  three  degrees  of  freedom.  By  fixing  temperature  and  amount  of  acid  in 
each  droplet,  one  degree  of  freedom  remains:  it  is  relative  humidity. 
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Figure  1 . Equilibrium  Diameter  of  Phosphoric  Acid  Droplets 
As  a Function  of  Relative  Humidity 

*Tarnove,  Theordore,  Munitions  Division,  Chemical  Systems  Laboratory,  Aberdeen  Proving  Ground, 
Maryland  21010,  February  1978,  private  communication. 
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Table  1 . Equilibrium  Diameter  of  Phosphoric  Acid  Droplets  As  a Function 
of  Relative  Humidity  (RH)  and  Condensation  Nuclei  Diameter  (dc) 
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each  droplet  is  made.  For  this  study,  the  mole  spectrum  of  H3PO4  varies  from  0.30  X 10“^  to 
0.14  X 10~12  moles.  It  will  sometimes  be  more  instructive  to  refer  to  a diameter  spectrum  for 
the  acid  nuclei,  denoted  by  dc;  and  for  this  study  dc  varies  from  0.30  to  2.5  pm.  These 
quantities  of  acid  will  be  considered  as  the  condensation  nuclei  for  further  vapor  diffusional 
growth.  With  reference  to  the  generation  of  phosphorus  smoke,  after  burning  the  elemental 
phosphorus  and  as  the  phosphorus  pentoxide  chemisorbs  water  vapor,  condensation  nuclei  are 
generated.  These  condensation  nuclei  consisting  of  phosphoric  acid  will  grow  by  vapor 
condensation  until  diffusional  equilibrium  with  the  environment  is  achieved.  The  “equilibrium 
states”  are  uniquely  specified  by  equations  (8)  and  (9). 

IV.  PHOSPHORIC  ACID  DROPLET  DIAMETER  AS  A FUNCTION  OF  RELATIVE 

HUMIDITY. 

At  this  point,  a relationship  between  saturation  ratio  over  the  surface  of  a droplet 
and  the  droplet  diameter  has  been  developed;  however,  as  stipulated  in  the  introduction,  the 
main  objective  is  to  obtain  a relationship  between  relative  humidity  and  droplet  diameter.  To 
achieve  the  stated  objective,  one  additional  statement  needs  to  be  furnished  and  that 
is:  diffusional  equilibrium  exists  between  the  surface  of  the  droplet  and  the  environment.  This 
statement  implies  equivalence  between  Se  and  Sj.  However,  Se  (the  environmental  saturation)  is 
synonymous  with  relative  humidity;  therefore,  through  transitive  association,  the  relationship 
between  Sj  and  d is  identical  to  the  relationship  between  RH  and  d,  and  the  desired  goal  has 
been  accomplished. 

With  this  final  identification,  the  predictive  relationships  between  relative  humidity 
and  phosphoric  acid  droplet  diameter  are  expressed  as, 


Ini// 


4a'Mj 

RTdp' 


12n2V, 


7rd3  - 6 n2(V2  + v ) 


1>  Xj  > .577 


(8a) 


ln<// 


4a  Mi 

7 ~ 

RTdp 


4.55  + 


27.3  n2V  \ 

7rd3  - 6 n2V2 


.577  > Xj  > 0 


(9  a) 


where  \ p is  the  relative  humidity  (RH).  After  specifying  the  constants  a',  Mj,  R,  T,  p',  Vj,  V2, 
and  the  moles  of  acid  (n2)  of  the  condensation  nucleus,  the  phosphoric  acid  droplet  diameter  is 
evaluated  as  a function  of  relative  humidity.  Therefore,  with  knowledge  of  the  condensation 
nuclei  size  distribution  of  an  oxidized  phosphorus  smoke,  the  phosphoric  acid  droplet  size 
distribution  can  be  predicted  from  the  environmental  relative  humidity. 

Scrutiny  of  table  1 reveals  the  following  salient  features  of  the  thermodynamic 
properties  of  H3PO4  droplets.  First,  the  relative  increase  in  droplet  diameter  is  independent  of 
the  initial  condensation  nucleus  diameter.  The  relative  diameter  increase  is  represented  by  the 
ratio  d^/dc,  where  d^  is  the  diameter  of  the  H3PO4  droplet  at  the  specified  relative  humidity 
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and  dc  is  the  diameter  of  the  particular  condensation  nucleus  instigating  growth.  From  the 
numerical  results  of  table  1,  calculations  for  d^/dc  at  various  \J/’s  reveal  that  the  ratio  is  only 
dependent  on  \ p and  not  on  dc.  Roughly,  the  relative  increase  in  droplet  diameter  takes  on  the 
values  1.5,  2,  and  3 for  relative  humidities  70%,  90%,  and  98%,  respectively , for  all  condensation 
nuclei  of  the  representative  phosphorus  smoke.  The  physical  reasoning  supporting  this  conclusion 
is  that  for  the  smoke  in  interest  the  solute  effect  dominates  the  Kelvin  effect.  Furthermore,  the 
solution  effect  acts  to  dilute  all  droplets  to  a uniform  extent.  Since  the  ratio,  d^/dc,  is  a 
parameter  of  dilution  and  all  droplets  achieve  the  same  degree  of  dilution  then  all  droplets 
experience  the  same  relative  increase  in  droplet  diameter,  irrespective  of  initial  condensation 
nucleus  size. 


Secondly,  the  change  in  droplet  diameter  concomitant  with  a change  in  relative 
humidity  is  a monotomically  increasing  function  of  relative  humidity  and  the  greatest  proportion  of 
growth  occurs  in  the  region  defined  by  > 80 %;  more  mathematically, 


9d/30  = H (\f)  t 


A cursory  mathematical  justification  for  the  “accelerated”  growth  of  droplets  with 
increasing  relative  humidity  can  be  developed  from  inspection  of  the  Kohler  expression.  For  high 
humidities,  the  phosphoric  acid  droplets  become  dilute  and  the  Kohler  expression  is  an  appropriate 
model  for  phosphorus  smoke.  Nevertheless,  recalling  that  Sj  is  equivalent  to  \jj  at  equilibrium,  the 
dependence  of  i//  on  d is  concluded  to  be  inversely  proportional.  Basic  differentiation  of  \p  with 
respect  to  d reveals  that  the  slope  of  this  expression,  3i///3d,  is  inversely  proportional  to  droplet 
diameter.  The  rate  of  change  of  droplet  diameter  with  respect  to  relative  humidity,  3d/3 \p,  obtained 
by  simple  inversion  of  3 i///3d  is  therefore  directly  proportional  to  droplet  diameter.  However, 
droplet  diameter  is  a monotomically  increasing  function  of  \p  and  3d/30  is  concomitantly  an 
increasing  function  of  relative  humidity,  \p. 

Finally,  it  is  noteworthy  to  realize  that  with  knowledge  of  the  droplet  diameter,  d,  and 
the  moles  of  acid,  n2,  in  the  droplet,  the  droplet  acid  concentration  can  be  calculated  once  the 
relative  humidity  is  specified  (see  appendix  C).  Information  on  the  droplet  acid  concentrations 
would  provide  some  information  on  the  optical  properties  of  the  phosphorus  smoke  since  acid 
concentration  is  tightly  linked  to  index  of  refraction. 

V.  APPROXIMATION  FORMULAE  FOR  PHOSPHORIC  ACID  DROPLET  DIAMETER  AS  A 

FUNCTION  OF  RELATIVE  HUMIDITY. 

As  formulated  in  Section  IV,  explicit  relationships  for  the  environmental  relative 
humidity  as  a function  of  droplet  diameter  and  number  of  moles  of  acid  (n2)  have  been  developed 
mathematically,  \p  = \p(d,  n2).  The  relationships  take  the  form: 


4a'Mi 
In  \p  = - 

RTp'd 
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1 + 


1 2 n2V , 

?rd3  - 6 n2(V2  + Vj) 


1 > Xj  > .577 
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4.55  + 


27.3  n2Vl 


(lib) 


4 o'M] 

Ini//  = In 

RTp'd 


7rd^  - 6 n2V2 


.577  > Xj  > 0 
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Clearly,  in  equations  (11a)  and  (lib)  the  droplet  diameter,  d,  assumes  the  role  of  the 
independent  variable  and  the  relative  humidity  \jj9  the  dependent  variable.  In  other  words,  by  first 
establishing  the  diameter,  the  relative  humidity  is  subsequently  known.  However,  in  reality,  the 
relative  humidity  assumes  some  value  independently  of  the  smoke  and  the  droplet  diameter  is 
established  once  \p  is  specified.  Therefore,  a functional  expression,  formatted  with  \jj  and  d as  the 
independent  and  dependent  variable  respectively,  is  physically  more  meaningful.  To  obtain  this 
expression,  equations  (1  la)  and  (lib),  transcendental  equations  in  diameter,  must  be 
mathematically  inverted. 

Unfortunately,  equations  (11a)  and  (lib)  cannot  be  inverted  analytically,  so  an 
approximate  explicit  relationship  for  the  droplet  diameter  must  suffice. 

As  concluded  in  Section  IV,  the  solution  effect  governs  the  thermodynamic  properties 
of  a phosphoric  acid  droplet  and  solely  determines  the  structure  of  the  relationship  between  droplet 
diameter,  d,  and  droplet  saturation  ratio,  S^.  In  phosphoric  acid  droplets  the  solution  effect  is 
mathematically  modelled  by 


Sd=2X]  - 1 


(12) 


where  Sd  and  X]  are  as  previously  defined. 

7rd^ 

Applying  the  substitutions,  Xj  = nj /(n j + n2)  and  =n2V2  + n]Vj,  to 

6 

equation  (1  2)  results  in 


Sd  = (7rd3  - 6 n2V2  - 6 n2  Vj  )/(7rd^  - 6 n2V,  + 6 n,V j ) 


(13) 


Equation  (13),  an  “equation  of  state”  for  the  properties  of  a phosphoric  acid  droplet,  provides  the 
defining  relationship  between  Sd  and  d. 


With  the  additional  substitution,  n2  = 7rpcdc'?/(6  M->),  where  pc  and  M2  are  the  density 
and  molecular  weight,  respectively,  of  phosphoric  acid,  a solution  for  d/dc  is  obtained  from 
equation  (13).  It  is 


d/dc  = (pc/M2) 


1/3 


V](l  +Sd)  + V2(l  -Sd)  1/(1  -Sd) 


1 1/3 


(14) 


For  the  phosphoric  acid-water  system,  the  following  data  are  relevant: 

pc  = 1.87  g/cc 
M2  = 98  g/cc 
Vj  = 1 8 cc/mole 
V2  = 52.4  cc/mole 
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Furthermore,  once  the  droplet  attains  equilibrium  with  the  environment,  is  equal  to  i p,  and 
equation  (14)  becomes 

r 1 1/3 

d/dc  = 1.1030  (1.  - 0.4886  i//)/(l.  - \p)  (15) 

Equation  ( 15)  is  an  approximate  explicit  relationship  for  the  phosphoric  acid  droplet  diameter  as  a 
function  of  relative  humidity  and  condensation  nucleus  diameter.  Comparison  of  exact  and 
approximate  diameters,  calculated  through  equations  (11)  and  (15),  respectively,  indicates,  at  most, 
a 2%  error.  The  extremely  close  agreement  corroborates  the  conclusion  that  the  solution  effect 
indeed  governs  the  properties  of  a phosphoric  acid  droplet  in  a phosphorus  smoke. 

Examination  of  equation  (15)  reveals  the  independence  of  the  relative  increase  in 
diameter  (d/dc)  on  the  condensation  nucleus  diameter  (dc).  Additionally,  d/dc  increases 
monotomically  with  \p  and  as  ^ -*•  1 the  increase  rises  sharply,  emphasizing  the  point  that  the 
greatest  proportion  of  growth  occurs  in  the  region  defined  by  \p  > 80%. 

Sensitivity  of  the  particle  size  distribution  function  on  relative  humidity  can  be 
investigated  with  knowledge  of  the  diameter  approximation  formula.  Consider  the  following 
illustration. 


As  elemental  phosphorus  is  burned,  particle  nucleation  occurs  giving  rise  to  a ln-normal 
distribution  for  the  condensation  nuclei.*  The  distribution  is  given  by 


where 


f(lndc)  = 


(2  7T)1/2 


exp 


o g 


[- 


(lndc-lndcg)-/2 


(16) 


f (In  dc)  = relative  number  of  condensation  nuclei  with 
diameters  between  In  d , and  In  dc  + A In  dc 

o.T  = geometric  standard  deviation 

o 

dc  = diameter  of  condensation  nucleus 
driy  = geometric  number  mean  diameter  of  distribution 


As  the  condensation  nuclei  grow,  the  original  distribution  changes  until  equilibrium  is 
reached  with  the  environmental  relative  humidity.  The  approximation  formula  contains  the 


*Test  chamber  experiments  on  oxidized  phosphorus  indicate  that  the  condensation  nuclei  are  characterized  by  a 
ln-normal  distribution. 
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information  that  determines  the  manner  in  which  the  distribution  varies.  By  simple  substitution  of 
equation  (15)  into  equation  (16),  the  new  size  distribution  of  droplet  diameters  is  given  by 


f(lnd) 


= £(2  ir)^2og 


-1 


exp 


-lnd  - ln(a(^/)dcg)  / 2 og2 


(17) 


where  a(i//)=  1.1030[(1.  - 0.4886i//)/(l . - i//)]1/3 

Equation  (17)  represents  the  size  distribution  function  of  a phosphorus  smoke  as  a 
function  of  relative  humidity  and  demonstrates  two  salient  points. 


1 . The  standard  deviation  (ag)  of  the  distribution  is  invariant  to  changes  in  relative 
humidity.  Therefore,  the  distribution  retains  its  structure  or  shape  through  condensational  growth. 


2.  The  central  tendency  or  geometric  mean  diameter  increases  from  dcg  to  a(<//)dcg 
and  results  in  a translation  of  the  size  distribution  toward  larger  droplet  diameters. 

VI.  KINETICS  OF  CONDENSATIONAL  GROWTH  OF  PHOSPHORIC  ACID  DRQP- 

LETS:  RELAXATION  TIMES. 


As  stated  earlier,  diffusional  equilibrium  is  established  when  the  environmental 
saturation  ratio  (relative  humidity)  is  equivalent  to  the  droplet  saturation  ratio.  Furthermore,  the 
linchpin  in  the  development  of  the  relationship  between  relative  humidity  and  droplet  diameter  is 
the  adoption  of  diffusional  equilibrium  between  the  environment  and  H3PO4  droplet.  To  test  the 
validity  of  this  assertion,  an  investigation  on  the  kinetics  of  condensational  growth  is  warranted. 
More  specifically,  development  of  a relaxation  time  will  provide  the  information  needed  to  answer 
the  question:  Is  diffusional  equilibrium  a reasonable  assumption  under  tactical  conditions? 

When  inequality  exists  between  the  saturation  ratios,  a gradient  in  vapor  concentration 
is  developed,  resulting  in  vapor  diffusion  toward  the  low  concentration  region.  Kinetic  growth  or 
decay  results  until  equality  of  the  saturation  ratios  is  once  again  established.  Before  launching  into  a 
quantitative  analysis  on  the  kinetics  of  condensational  growth,  it  is  instructive  to  illustrate  some 
qualitative  features  of  this  nonequilibrium  phenomenon.  More  specifically,  with  information  on  the 
properties  of  a droplet  and  the  environmental  relative  humidity,  conclusions  concerning  growth  or 
decay  may  be  ascertained. 

With  reference  to  figure  2,  a representative  equilibrium  curve  is  diagrammed  with  the 
impressed  external  condition  of  an  unsaturated  environment,  i.e. , Se  < 1 . For  convenience,  consider 
a continuous  distribution  of  droplets  along  the  curve,  each  with  its  characteristic  saturation  ratio 
defined  by  the  curve.  Then  one,  and  only  one,  droplet  along  this  curve  will  be  in  equilibrium  with 
the  environmental  relative  humidity  defined  by  the  intersection  of  Se  and  Sj.  This  intersection 
represents  the  condition  for  equality  of  the  saturation  ratios  which  is  equivalent  to  diffusional 
equilibrium.  For  those  droplets  distributed  along  the  lower  portion  of  the  curve,  below  Se,  the 
characteristic  droplet  saturation  ratios  are  less  than  the  environmental  saturation  ratio.  Vapor 
diffusion  toward  the  droplets  occurs,  resulting  in  condensational  growth  until  their  diameters 
correspond  to  the  equilibrium  diameter  for  that  environment.  For  those  droplets  distributed  along 
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Figure  2.  Kinetics  of  Approach  to  Equilibrium  by  an 
Aerosol  in  an  Unsaturated  Environment 


20 


the  upper  portion  of  the  curve,  Sd  > Se,  evaporation  will  result  until  they  arrive  at  the  equilibrium 
diameter,  aj . It  is  a simple  matter  to  envision  a family  of  curves,  each  possessing  an  “equilibrium” 
diameter,  thereby  generating  a collection  of  equilibrium  diameters  for  the  given  environment. 

The  same  analysis  may  be  applied  to  figure  3 where  the  external  constraint  is  a 
supersaturated  environment.  Now  two  equilibrium  diameters  are  present  (two  intersections):  <Xj 
and  «2-  However,  these  equilibrium  diameters  are  distinct  not  only  in  size  but  also  in  nature.  The 
smaller  diameter,  aj , is  designated  a stable  equilibrium  point,  whereas  o :2  is  referred  to  as  an 
unstable  equilibrium  diameter.4  The  basis  for  these  distinctions  may  be  explained  as 
follows:  Consider  a droplet  with  diameter  ctj , if  displaced  upward  along  the  curve,  holding  Se 
constant,  an  evaporation  process  shall  result  until  equilibrium  is  once  established  at  <Xj.  Conversely, 
if  displaced  downward,  then  condensation  occurs  until  equilibrium  is  restored.  Therefore,  with  any 
displacement  from  equilibrium  diameter  a restoring  force  is  established  in  the  form  of 
differential  saturation  ratios.  The  same  analysis  reveals  that  for  any  displacement  from  a2,  a driving 
force  away  from  equilibrium  is  established,  thereby  arriving  at  the  nomenclature  of  an  unstable 
equilibrium  diameter.  Of  special  interest  is  the  phenomenon  of  free  growth  resulting  from  a 
displacement  away  from  the  equilibrium  diameter  a-,  toward  larger  diameters.  When  this  occurs, 
holding  Se  constant,  the  droplet  may  freely  grow  by  condensation,  resulting  in  unbounded  growth 
toward  “giant”  diameters.  For  the  representative  phosphorus  smoke  the  supersaturation  ratio  is 
only  .12%,  an  extremely  small  value  for  the  initiation  of  free  growth. 

To  further  the  understanding  of  condensational  kinetics,  a quantitative  analysis  is  in 
order.  Information  on  the  time  for  approach  to  equilibrium  is  of  considerable  importance  and 
illuminates  the  problem  of  diffusional  equilibrium.  The  parameter  which  characterizes  this  time  is 
designated  the  relaxation  time,  r,  and  may  be  introduced  through  the  modified  diffusional  growth 
law  equation  (18). 


(a  + £)  da/dt  = p Deff  (se  “ S(a)/P£ 
eq 


(18) 


where 


a = instantaneous  radius  of  the  droplet 

£ and  De«-  = surface  effect  parameters  introduced  to  generalize  the  growth  law  for 
the  transition  region  between  molecular  and  continuum  phenomenon 

p — equilibrium  vapor  density  of  the  environment  at  that  temperature 
eq 

pg  = density  of  the  diffusing  vapor 
t = time 
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Figure  3.  Kinetics  of  Approach  to  Equilibrium  by  an 
Aerosol  in  a Super  Saturated  Environment 


The  modified  growth  law  developed  by  Carstens  et  al.4  represents  a balance  between 
vapor  diffusion  and  heat  conduction  enveloped  within  a quasi-steady  state  analysis.  Expressions  for 
the  surface  effect  parameters  have  been  developed  by  Carstens  involving  the  condensation 
coefficient  a and  the  evaporation  coefficient/?.  For  this  study,  the  values  a = 1.0  and  (3=  0.5  have 
been  employed. 

The  solution  to  the  first-order  differential  equation,  a(t),  is  in  general  difficult  to 
achieve  analytically.  However,  for  the  constraint  growth  toward  stable  equilibrium,  an  approximate 
analytical  solution  may  be  displayed. 

Mathematical  reforming  of  equation  (18)  results  in 


where 


a(a'  + g)da' 

Se-S(a') 


(19) 


aQ  = initial  radius 

t (a)  = time  at  which  the  droplet  achieves  final  radius  a 


With  substitution  of  equations  (8)  and  (9)  for  S(a)  in  equation  (19),  the  solution  of  the  problem 
involves  the  integration  of  a complicated  function  of  a'.  An  exact  integration  is  tedious,  and 
approximate  integration  techniques  are  commonly  employed.3  (See  appendix  A.)  With  the 
aforementioned  integration  techniques,  a solution  a(t)  is  obtained  in  the  form: 


(a 


-aj)  = (a0-  cq)e  t,Tal 


(20) 


where 


cq  - equilibrium  radius  defined  by  Se 

tQ!j  = relaxation  time,  characteristic  time  for  approach 

to  equilibrium  as  growth  occurs  from  aQ  to  a, . NOTE: 

t -*  °°  a -*■  cq 

and 

t - 0 a = a0  as  expected. 

The  exponential  dependence  on  time  for  a(t)  is  a result  of  the  basic  form  of  the  growth 
law  and  the  condition  of  growth  toward  stable  equilibrium.  To  gain  further  insight  into  this 
solution,  a graphical  representation  has  been  developed  (figure  4).  We  see  a droplet  with  initial 
radius  a0  (aQ  is  the  radius  of  the  condensation  nucleus)  growing  by  condensation  toward  the 
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Figure  4.  Approach  to  Stable  Equilibrium  for  an  Unsaturated  Environment 


equilibnum  radius  . The  instantaneous  radius  a(t)  is  defined  by  equation  (20)  involving  the 
relaxation  time,  ra^,  which  represents  the  time  for  approach  to  . The  analytical  form  of  Ta  is  in 

appendix  A and  takes  the  form 


a 


= K(«j  + fi)  (8  7raj4  - 21 1.8  n2a1 )/  ^8  tt(  1 - \p)  (oij  - a2)  (ai  _ “3)  (“1  - <*4)]  (21) 


The  numerical  value  of  Ta  rests  upon  the  values  K,  £,  n2,  \p,  , a2,  (*3,  and  a4>  The  parameters,  K 

and  £,  can  be  considered  a priori  constants  independent  of  variations  in  \ p and  d.  The  amount  of 
acid,  n2,  in  the  droplet  of  interest  must  be  independently  specified  thereby  acting  as  a constraint  on 
growth.  Once  \p  is  specified,  the  equilibrium  diameter  a 1 along  with  a2  and  a 4 are  established. 
The  variables  a.^  a^,  ^3?  and  represent  solutions  to  the  quartic  expression  in  radius  (see 
appendix  A,  equation  (3A))  and  are  derivable  once  \Jj  is  established.  Consequently,  as  \p  and  n2  are 
named,  all  other  variables  in  equation  (21)  become  known  and  calculation  of  is  completed. 

Computations  for  ra  are  performed  for  values  of  \p  (RH)  ranging  from  10%  to  98%  and  n2  ranging 
from  approximately  0.30  E-15  to  0.74  E-12  moles  of  acid,  and  are  displayed  in  table  2. 


Fortunately,  the  growth  toward  stable  equilibrium  typifies  the  growth  of  phosphorus 
smoke  under  tactical  conditions  and  therefore  equation  (20)  should  represent  an  adequate  model 
for  the  initial  formation  of  a phosphorus  smoke  cloud.  Specifically,  aQ  represents  the  radius  of  the 
condensation  nucleus  introduced  earlier,  and  with  an  analytical  expression  for  ra  , a tabular 

presentation  for  relaxation  time  versus  relative  humidity  is  displayed  in  table  2. 


As  seen  from  table  2,  numerical  results  reveal  extremely  small  relaxation  times  for 
environments  characterized  by  relative  humidity  ranging  from  10%  to  98%.  Orders  of  magnitude  for 
the  relaxation  times  range  from  10 to  0.90  seconds. 


In  other  words,  for  the  case  of  unsaturated  environments,  acid  droplets  attain 
equilibrium  with  the  environment  almost  instantaneously , and  the  assertion  of  diffusional 

equilibrium  under  tactical  conditions  is  valid.  * As  a result , acid  droplets  can  be  said  to  “track"  with 
the  environment , a condition  conducive  to  a strong  relationship  between  droplet  properties  and 

en  v iro  nmental  prop  erti  es . 


As  developed  in  appendix  B,  the  relaxation  time  r is  inversely  proportional  to  the  slope 
of  the  saturation  ratio  curve,  evaluated  at  the  point  of  equilibrium.  Mathematically, 


where 


(22) 


S = S(a) 


^Assuming  the  time  of  combustion,  time  of  chemisorption,  and  time  of  hydration  is  negligible. 
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Table  2.  Relaxation  Times  for  Approach  to  Equilibrium  As  a Function 
of  Relative  Humidity  and  Moles  of  H3PO4  (N) 
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With  this  identification  for  t,  it  is  realized  that  as  the  slope  of  the  curve  becomes  large  (steep),  the 
relaxation  time  becomes  small;  and  with  inspection  of  figure  1,  one  concludes  that  characteristic 
relaxation  times  for  unsaturated  environments  (RH  < 100)  are  of  extremely  small  magnitude. 
Another  qualitative  feature  of  table  2 is  the  strong  correlation  between  moles  of  solute  and 
relaxation  time,  t.  Specifically,  the  larger  the  condensation  nuclei  (moles  of  solute)  the  greater  the 
relaxation  time,  r,  indicating  a longer  period  for  growth  toward  stable  equilibrium.  This 
circumstance  is  evidenced  in  equation  (19)  where  the  time  t(a)  is  inversely  proportional  to  the 
difference  Se  - S(a).  The  difference  is  greater  for  droplets  containing  larger  moles  of  solute, 
resulting  in  longer  periods  of  time  for  approach  to  equilibrium.  An  alternate  perspective  is  obtained 
from  the  dilution  effect.  In  equilibrium,  all  droplets  tend  to  have  the  same  acid  concentration.  For 
larger  moles  of  solute,  a greater  addition  of  water  is  necessary  in  order  to  achieve  the  same  degree  of 
dilution  and  therefore  the  greater  relaxation  times. 


With  a valid  relationship  between  relative  humidity  and  equilibrium  drop  size,  yield 
factors  for  a phosphorus  smoke  may  be  ascertained. 

VII.  YIELD  FACTOR  VERSUS  RELATIVE  HUMIDITY  FOR  PHOSPHORUS  SMOKE. 

Yield  factor  may  be  defined  as  the  ratio  of  final  mass  to  initial  mass.  In  the  case  of 
phosphorus  smoke,  the  initial  mass  corresponds  to  the  amount  of  red  or  white  phosphorus  burned, 
and  the  final  mass  is  determined  by  the  mass  content  of  the  smoke  cloud.  The  mass  of  the  smoke 
cloud  is  the  result  of  four  basic  reactions:  combustion,  chemisorption,  hydration,  and 
condensation.9 

Each  reaction  represents  an  addition  of  mass  to  the  initial  phosphorus  burned. 
Combustion  results  in  an  addition  of  oxygen  whereas  chemisorption,  hydration,  and 
condensation  add  water  to  the  system.  The  estimated  size  of  the  condensation  nuclei  accounts  for 
the  addition  of  oxygen  and  water  due  to  combustion,  chemisorption,  and  hydration.  The 
specification  of  drop  size  through  the  developed  model,  with  knowledge  of  condensation  nuclei, 
theoretically  accounts  for  water  addition  through  condensation  and,  therefore,  the  determination  of 
the  mass  of  the  smoke  cloud  is  uniquely  determined  by  the  drop  sizes  defined  by  the  model. 
Obviously,  drop  sizes  are  relative  humidity-dependent  and  consequently  yield  factors  will  be  relative 
humidity-dependent. 

To  initiate  the  quantitative  analysis  for  yield  factors,  consider  the  following  three-step 
reaction  accounting  for  combustion,  hydration,  and  condensation. 


d , rr\  combustion  chemisorption 

4(s)  5%,  > P4°>°(V),(S)  +6H2°«  »4H3P°4 


(2) 

(condensation 

nuclei) 


(23) 


+ nH^O 


condensation  + hydration 


»h3po4 


(eq) 


(drop  sizes  determined  by  model) 
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Therefore,  for  each  P4  complex  burned,  a (4H3P04)  complex  is  generated  representing  the 
condensation  nucleus  for  further  growth.  From  a mass  bookkeeping  point  of  view,  generation  of  the 
smoke  cloud  is  equivalent  to  generation  of  the  condensation  nuclei,  nature  “picking  up  the  mass 
tab”  for  condensation.  In  order  to  generate  the  ith  condensation  nucleus,  consisting  of  4H3P04 
complexes,  a specific  amount  of  P4  must  be  burned.  The  following  algorithm  is  used:  let  the  mass 
of  the  ith  condensation  nucleus  be  represented  by  Mj,  then  the  number  of  4H3P04  complexes 
contained  in  the  ith  condensation  nucleus,  denote  it  Nj,  is  given  by 


M(4H3P04) 


(24) 


Since  for  each  4H3P04  complex  one  P4  complex  must  be  burned,  then  Nj  is  identically  the  number 
of  P4  complexes  needed  to  generate  the  ith  condensation  nucleus.  This  establishes  the  initial  mass 
burned  needed  to  generate  the  ith  condensation  nucleus. 

Therefore,  the  yield  factor  due  to  the  ith  droplet  (CN  + growth)  is  given  by: 


Y:  = — a3 


N:Mp 
1 1 4 


where 


(25) 


peff  = (%  H3P04  + % H20 )/(%  H3P04/pH3P04  + % H20/pH20) 


where 


aj  = radius  of  the  ith  droplet  at  the  specified  relative  humidity 
% H3P04  and  % H20  = percent  acid  and  water  in  the  droplet  at  that  relative  humidity 

Finally,  with  knowledge  of  fj  representing  the  relative  frequency  of  the  ith 
condensation  nucleus,*  by  taking  a weighted  average  over  all  condensation  nuclei,  the  total  yield 
factor  as  a function  of  relative  humidity  may  be  written  as 

YF(RH)  = Sf^rra]  (peff/3NiMP4j  (26) 

The  numerical  results  are  displayed  graphically  in  figure  5 and  table  3.  Yield  factors 
range  from  3.8  to  16  for  RH’s  ranging  from  10%  to  90%.  The  greatest  proportion  of  growth  occurs 
in  the  region  defined  by  RH  > 80%.  The  augmentation  of  mass  is  considerable  and  indicates  the 
importance  of  condensation  in  the  determination  of  yield  factors  for  hygroscopic  smokes  such  as 
phosphorus  smoke. 
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Figure  5.  Yield  Factors  of  a Phosphorus  Smoke  As 
a Function  of  Relative  Humidity 
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Table  3.  Yield  Factors  of  a Phosphorus  Smoke  As 
a Function  of  Relative  Humidity 


YF 

3.82 

4.42 

4.91 

5.04 

5.36 

5.97 

7.13 

9.33 

16.29 

RH 

10% 

20% 

35.8% 

40% 

50% 

60% 

70% 

80% 

90% 

VIII.  CONCLUSIONS. 

It  has  been  concluded  that: 

1.  With  knowledge  of  the  condensation  nuclei  size  distribution  of  an  oxidized 
phosphorus,  the  droplet  size  distribution  of  a phosphorus  smoke  can  be  confidently  predicted. 

2.  The  relative  increase  in  droplet  diameter  is  solely  dependent  on  relative  humidity 
and  is  uniform  for  all  acid  droplets  in  a phosphorus  smoke.  The  relative  increase  in  diameter, 
referenced  to  the  condensation  nucleus  diameter,  is  1.5,  2,  and  3 for  relative  humidities  70%,  90%, 
and  98%. 


3.  The  change  in  droplet  diameter  concomitant  with  a change  in  relative  humidity 
is  a monotomically  increasing  function  of  relative  humidity  and  the  greatest  proportion  of  growth 
occurs  in  the  region  defined  by  RH  >80%. 

4.  Droplet  acid  concentrations  of  a phosphorus  smoke  can  be  effectively 
ascertained  through  the  size  prediction  model,  thereby  providing  information  on  the  optical 
properties  of  a phosphorus  smoke. 

5.  The  diameter  of  a phosphoric  acid  droplet  can  be  explicitly  expressed  as  a 
function  of  relative  humidity  and  condensation  nucleus  diameter  with  sufficient  accuracy. 

6.  Relaxation  times  for  approach  to  equilibrium  under  condensational  growth  are 
generally  less  than  1 second  for  a phosphorus  smoke;  consequently,  the  phosphoric  acid  droplets 
should  "track”  with  the  environment  under  tactical  conditions. 

7.  Theoretical  yield  factors  are  nonlinearily  dependent  on  relative  humidity  and 
vary  from  3,89  to  16.29  at  relative  humidities  10%  and  90%,  respectively,  for  a phosphorus  smoke. 
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GLOSSARY 


o 

a',b' 

d 


Jcg 

eff 


D, 

M] 

M; 


SO 


R 

RH,  * 


T 

t 

V1 

v2 

X, 

Xleff 

YF 

o' 

o 

z 

p 

pZ 

p(°°) 

H eq 
Peff 


Solution  droplet  radius 
Initial  radius  of  solution  droplet 
Empirical  functions  of  relative  humidity 
Solution  droplet  diameter 
Condensation  nucleus  diameter 

Geometric  mean  diameter  of  condensation  nuclei  size  distribution 
Effective  diffusion  coefficient;  a surface  effects  parameter 
Molecular  weight  of  solvent  (water) 

Mass  of  ith  condensation  nucleus 
Moles  of  solvent  (water) 

Moles  of  solute  (H3PO4) 

Actual  environmental  vapor  pressure  (water) 

Equilibrium  vapor  pressure 
Vapor  pressure  over  droplet 
Universal  gas  constant 
Relative  humidity 
Droplet  saturation  ratio 
Environmental  saturation  ratio 
Temperature 
Time 

Mole  volume  of  solvent 
Mole  volume  of  solute 
Mole  fraction  of  solvent 
Effective  mole  fraction  of  solvent 
Yield  factors 

Surface  tension  of  droplet 

Standard  deviation  of  size  distribution 

Equilibrium  diameter 

Surface  effects  parameter 

Density  of  bulk  solution 

Density  of  diffusing  vapor 

Equilibrium  vapor  density  of  environment 

Effective  density  of  solution  droplet 
Relaxation  time 
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APPENDIX  A 


GROWTH  INTEGRAL 


Approximate  integration  of  growth  integral  for  approach  to  stable  equilibrium; 
relaxation  time  r. 

a 


With 


t(a)  = 


/“  (a+g)da 
„ Se  - S(a) 


(1A) 


eq  “o 

substitution  for  S(a)  given  by  equation  6 results  in  the  following  integrand  involving  the  variable  a 
only: 


t(a)  = K 


I 


(a'  + C)da' 


where 


0 Se-e 


.23  X 10-6/d  - £n  (1  + 


216n2 


7raJ  - 427.8  n2  /_ 


K = 


p (oo)  Deff 
eq 


Now,  since  .23  X 10  ^/d«l  for  diameters  of  interest,  the  approximation  ex  = 1 + x may  be 
applied  to  S(a)  resulting  the  form, 

r 1 


S(a)  = (1  + .23  X lO^/d) 


1 + • 


2l6n2 


7rd3  -427.8n 


2 J 


After  standardizing  S(a)  into  a proper  fraction,  the  time  integrand  becomes, 

a 

(a'  + £)da 


t(a)  = K 


I 


(a + .11  N 10-6)  (87ra3  -319.8n2) 


ao 


36.4 7ra4  - 963.69n2a 
and  with  the  additional  standardization  of  the  time  integrand  we  arrive  at; 


t(a)  = K 


/ 


(a' + C)  (87ra4  - 21 1.8n2a)da' 


87T (S  — l)a4  - .88  X lO"6™3  + (427.8n2  - 21 1.8n2Sa  + 47.06  X 10"6n2 


(2A) 


To  evaluate  this  complex  integrand,  follow  the  given  procedure: 

By  factoring  87r(s  - 1)  from  the  denominator,  a quartic  expression  (4th  order  polynomial  in  a) 
is  realized  where; 


35 


t(a)  = 


K 


(a'  + fi)  (87ra4-211.8n2a)da 


8tt(S 


-f 

+*  a 


a4  — 


.88  X 1Q-6tt  3 t / 427.8  n2  - 21 1 .8n2s\  47.06  X 10~6 


8tt(s-1) 


n2 


8tt(S  - 1) 


8tt(S-1) 


Now  the  quartic  expression  a4  + ba3  + ca  + d may  be  rewritten  as  the  product  of  its  factored  roots, 
i.e., 


a4  + ba3 


+ ca  + d = (a  - oq ) (a  - a2)  (a  - oq)  (a  - a4) 


where  oq , «2,  «3,  and  a4  are  obtained  by  solving  the  equation  a4  + ba3  + ca  +d  = 0.  (This  is  easily 
handled  by  most  pocket  calculators.)  Therefore,  the  time  integrand  becomes: 

• a 


t(a)  = 


K 


8tt(S-1) 


(a  +•£)  (87ra4-211.8n2a)da 
(a-aq)  (a-a2)  (a -a3)  (a-a4) 


(3  A) 


where  oq , a2,  a3,  and  «4  satisfy  the  polynomial  equation, 

87r(S- l)a4  - .88  X 10-b7ra3  + (427. 8n2  - 21 1.8n2S)a  + 47.06  X 10_6n2  = 0. 


Now  t(a)  represents  the  time  at  which  the  droplet,  experiencing  condensational  growth,  achieves  its 
final  radius  a.  For  the  problem  of  approach  to  stable  equilibrium,  the  final  radius  is  the  radius  of 
equilibrium  and  is  given  by  the  condition  that  Sg  = S(a).  This  relationship  is  satisfied  by  the  values 
of  a such  that  a = oq , a = «2,  a = a3,  and  a = «4.  Now  only  one  of  these  four  solutions  is  positive 
and  real,  call  it  oq,  and  it  represents  the  equilibrium  radius.  Thus  the  upper  limit  of  integration  in 
equation  (3 A)  is  identical  to  the  final  equilibrium  radius  aq . However,  when  a = oq  , the 
denominator  of  the  integrand  in  equation  (3A)  is  zero  and  the  integrand  is  consequently  infinite.  In 
other  words,  due  to  the  singular  nature  of  l/(a-oq),  the  magnitude  of  the  integrand  is  sharply  peaked 
in  the  neighborhood  of  a = oq , as  seen  in  the  figure.  From  a time  viewpoint,  this  implies  that  for 
growth  toward  stable  equilibrium  the  greatest  portion  of  time  is  accounted  for  in  the  neighborhood 
of  equilibrium.  This  statement  is  confirmed  by  the  modified  growth  law  which  states  that  the 
change  in  radius  is  greatest  when  the  difference  in  saturation  ratio  is  greatest.  Near  equilibrium  we 
have  a minimum  in  this  difference  and,  therefore,  a minimum  in  the  rate  of  change  in  radius. 
Therefore,  a long  time  is  spent  in  condensational  growth  near  equilibrium.  Thus,  to  a good 
approximation,  all  terms  in  the  integrand  are  mathematically  constant  as  compared  to  the  highly 
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singular  term  l/(a  - c^).  Therefore,  a mean  value  approach  is  in  order  and  all  terms  may  be  brought 
outside  the  integral  evaluated  at  aj  except  for  the  term  contributing  to  the  largeness  of  the 
integrand  (l  /(a  - c^)) . Therefore,  we  may  mathematically  write: 


t (a)  ~ 


K (oq  +£)  (87ra!4  -211.8n2o:i) 


87r(S- 1)  (aj -a2)  (ai -0;3)  (ai -a4)  I (a-ai) 


/ 


a da' 


Now  S - 1 < O for  unsaturated  environments  and  we  may  rewrite  t(a)  as, 

✓ ^ , a-al 

t(a)  — “ ro;jln  ■ 

where 


ao-al 


'a 


1 


K(aj+£)  (87raj4  - 21 1.8n2aj) 

87T ( IS-  1 1)  (aj  -a2)  (aj  -a3)  (aj  -«4) 


(4A) 


(5A) 


Now  equation  (4A)  may  be  reformed  into  its  exponential  equivalent  giving  us  the  growth  solution 
to  stable  equilibrium  characterized  by  the  relaxation  time  ra  given  by  equation  (5A),  where 


(a-aj)  = (a0"al)e  t/tai 
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APPENDIX  B 


ANALYTICAL  EXPRESSION  FOR  RELAXATION  TIME 


With 


(a'  + £)daf 
Se-S(a) 


and  the  understanding  that  this  integrand  is  dominated  in  the  vicinity  of  a = aj , then  we  may 
approximate  this  integrand  in  the  following  manner: 


Since  this  integrand  is  dominated  near  a = o(j,  we  can  expand  S(a)  in  a Taylor  series 
expansion  about  retaining  only  the  first  term, 


Therefore 


and  therefore 


But 


S(a)  = S(a1)  + /-—  J (a-a1)  + l 
\aa/a  = a1 


t(a)  = K 


/ 


(a'  + £)da' 


ao  ^e  ^(ai) 


_as 

3a 


a = a 


f a 

1 


S^)  = Se  t(a)  = K 

Jo  V 3a  ) _ 

N 7a  = aj 

Again  using  the  mean  value  approach  of  appendix  B we  obtain 


(a  - otj ) 


+ £)da' 


(a  — ocj ) 


T 


“1 


K(«!  +£) 


and,  therefore,  the  relaxation  time  is  inversely  proportional  to  the  slope  of  the  saturation  curve 
evaluated  at  the  radius  of  equilibrium. 
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APPENDIX  C 


ACID  CONCENTRATIONS  OF  PHOSPHORIC  ACID  DROPLETS 

AS  A FUNCTION  OF  RELATIVE  HUMIDITY 


Relative  humidity 

H3PO4  Acid 

% 

% 

98 

6 

90 

23 

80 

38 

70 

49 

60 

58 

50 

65 

40 

70 

36 

72 

20 

78 

10 

82 

0 

100 
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DISTRIBUTION  LIST  6 


Names 

CHEMICAL  SYSTEMS  LABORATORY 

SAFETY  OFFICE 
Attn:  DRDAR-CLF 
PLANNING  & PROGRAMS  OFFICE 
Attn:  DRDAR-CLR-L 
Attn:  DRDAR-CLR-T 

AUTHOR’S  COPIES:  Research  Division 

BIOMEDICAL  LABORATORY 
Attn:  DRDAR-CLL-B 
Attn:  DRDAR-CLL-M 
Attn:  DRDAR-CLL-MC 
Attn:  DRDAR-CLL-V 

DEVELOPMENTAL  SUPPORT  DIVISION 
Attn:  DRDAR-CLJ-R 
Attn:  DRDAR-CU-L 
Attn:  DRDAR-CLJ-M 
Attn:  DRDAR-CLJ-P 

MUNITIONS  DIVISION 
Attn:  DRDAR-CLN 

RESEARCH  DIVISION 
Attn:  DRDAR-CLB 
Attn:  DRDAR-CLB-B 
Attn:  DRDAR-CLB-C 
Attn:  DRDAR-CLB-P 
Attn:  DRDAR-CLB-R 
Attn:  DRDAR-CLB-T 
Attn:  D RD  AR-CLB-TE 

SYSTEMS  ASSESSMENTS  OFFICE 
Attn:  DRDAR-CLY-A 
Attn:  DRDAR-CLY-R 

DEPARTMENT  OF  DEFENSE 

Administrator 

Defense  Documentation  Center 

Attn:  Document  Processing  Division  (DDC-DD) 
Cameron  Station 
Alexandria,  VA  22314 

Director 

Defense  Intelligence  Agency 
Attn:  DB-4G1 
Washington,  DC  20301 

DEPARTMENT  OF  THE  ARMY 

HQDA  (DAMO-SSC) 

WASH  DC  20310 

Director 

Defense  Civil  Preparedness  Agency 
Attn:  PO(DC) 

Washington,  DC  20301 


Names  Copies 

CINCUSAREUR 

Attn:  AEAGC-RSI  1 

APO  New  York  09403 

Deputy  Chief  of  Staff  for  Research, 

Development  & Acquisition 

Attn:  DAMA-CSM-CM  1 

Attn:  DAMA-ARZ-D  1 

Washington,  DC  20310 

US  Army  Research  and  Standardization 

Group  (Europe)  1 

Attn:  Chief,  Chemistry  Branch 
Box  65,  FPO  New  York  09510 

HQDA  (DAMI-FIT)  1 

WASH,  DC  20310 

Commander 

HQ  US  Army  Medical  Command,  Europe 

Attn:  AEMPM  1 

APO  New  York  09403 


US  ARMY  HEALTH  SERVICE  COMMAND 

Superintendent 
Academy  of  Health  Sciences 
US  Army 

Attn:  HSA-CDC  1 

Attn:  HSA-IHE  1 

Fort  Sam  Houston,  TX  78234 

US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Attn:  DRCLDC  1 

Attn:  DRCSF-P  1 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 

Office  of  the  Project  Manager  for  Chemical  Demilitarization 
and  Installation  Restoration 

Attn:  DRCPM-DR-T  2 

Aberdeen  Proving  Ground,  MD  21010 

Project  Manager  Smoke /Obscurance 

Attn:  DRCPM-SMK-M  1 

Aberdeen  Proving  Ground,  MD  21005 

Human  Engineering  Laboratory  HFE  Detachment 

Attn:  DRXHE-EA  1 

Building  E3220 

Aberdeen  Proving  Ground,  MD  21010 
Commander 

US  Army  Foreign  Science  & Technology  Center 

Attn:  DRXST-MT2  1 

220  Seventh  St.,  NE 
Charlottesville,  VA  22901 


Copies 

1 

6 

1 

25 

1 

1 

1 

1 

3 

3 

1 

1 

5 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

1 

1 

1 
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DISTRIBUTION  LIST  6 (Contd) 


Names  Copies 

Director 

DARCOM  Field  Safety  Activity 

Attn:  DRXOS-C  1 

Charlestown,  IN  47111 

Commander 

US  Army  Armament  Materiel  Readiness  Command 

Attn:  DRSAR-ASN  1 

Attn:  DRSAR-SA  1 

Rock  Island,  IL  61299 

US  ARMY  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  COMMAND 

Commander 

US  Army  Armament  Research  and  Development  Command 


Attn:  DRDAR-AC  1 

Attn:  DRDAR-LCA  1 

Attn:  DRDAR-LCU  1 

Attn:  DRDAR-SCA-C  1 

Attn:  DRDAR-SCP-A  1 

Attn:  DRDAR-SCW  1 

Attn:  DRDAR-SER  1 

Attn:  DRDAR-TSS  2 

Dover,  NJ  07801 

Director 

Ballistic  Research  Laboratory 

Attn:  DRDAR-TSB-S  2 

Building  305 


Aberdeen  Proving  Ground,  MD  21005 

CDR,  APG 
USA  ARRADCOM 

Attn:  DRDAR-GCL  1 

Aberdeen  Proving  Ground,  MD  21010 

US  ARMY  ARMAMENT  MATERIEL  READINESS  COMMAND 
Commander 

US  Army  Armament  Materiel  Readiness  Command 


Attn:  DRSAR-IMB-C  1 

Attn:  DRSAR-PDM  1 

Attn:  DRSAR-SA  1 

Attn:  DRSAR-SF  1 

Rock  Island,  1L  61299 

CDR,  APG 
USA  ARRCOM 

Attn:  DRSAR-MAS-C  1 

Attn:  SARTE  1 


Aberdeen  Proving  Ground,  MD  21010 
Commander 

US  Army  Dugway  Proving  Ground 

Attn:  Technical  Library,  Docu  Sect  1 

Dugway,  UT  84022 

Commander 

Rocky  Mountain  Arsenal 

Attn:  SARRM-OA  1 

Commerce  City,  CO  80022 


Names  Copies 

Commander 
Pine  Bluff  Arsenal 

Attn:  SARPB-ETA  1 

Pine  Bluff,  AR  71611 

US  ARMY  TRAINING  & DOCTRINE  COMMAND 

Commandant 

US  Army  Infantry  School 

Attn:  NBC  Division  1 

Fort  Benning,  GA  31905 


Commandant 

US  Army  Missile  & Munitions  Center  & School 

Attn:  ATSK-CD-MD  1 

Attn:  ATSK-DT-MU-EOD  1 

Redstone  Arsenal,  AL  35809 

Commandant 

US  Army  Military  Police  School/Training  Center 

Attn:  ATZN-CDM  1 

Attn:  ATZN-TDP-C  4 

Fort  McClellan,  AL  36205 

Commander 

US  Army  Infantry  Center 

Attn:  ATSH-CD-MS-C  1 

Fort  Benning,  GA  31905 

Commandant 

US  Army  Ordnance  & Chemical  Center  & School 

Attn:  ATSL-CL-CD  1 

Aberdeen  Proving  Ground,  MD  21005 


US  ARMY  TEST  & EVALUATION  COMMAND 
Commander 

US  Army  Test  & Evaluation  Command 

Attn:  DRSTE-FA  1 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Cold  Regions  Test  Center 

Attn:  STECR-TD  1 

APO  Seattle,  WA  98733 

DEPARTMENT  01  THE  NAVY 

Commander 

Naval  Explosive  Ordnance  Disposal  Facility 

Attn:  Army  Chemical  Officer,  Code  604  1 

Indian  Head,  MD  20640 

Commander 

Nuclear  Weapons  Training  Group,  Atlantic 
Naval  Air  Station 

Attn:  Code  21  1 

Norfolk,  VA  23511 

Chief,  Bureau  of  Medicine  & Surgery  1 

Department  of  the  Navy 
Washington,  DC  20372 
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DISTRIBUTION  LIST  6 (Contd) 


Names  Copies 

Commander 

Naval  Weapons  Center 

Attn:  A.  B.  Galloway /Code  3171  1 

China  Lake,  CA  93555 

Commanding  Officer 

Naval  Weapons  Support  Center 

Attn:  Code  5042/Dr.  B.  E.  Douda  1 

Crane,  IN  47522 

US  MARINE  CORPS 

Director,  Development  Center 

Marine  Corps  Development  & Education  Command 

Attn:  Fire  Power  Division  1 

Quantico,  VA  22134 

DEPARTMENT  OF  THE  AIR  FORCE 

HQ  Foreign  Technology  Division  (AFSC) 

Attn:  PDRR  1 

Wright-Patterson  AFB,  OH  45433 

Commander 

Aeronautical  Systems  Division 

Attn:  ASD/AELD  1 

Wright-Patterson  AFB,  OH  45433 

HQ,  USAF/SGPR  1 

Forrestal  Bldg 
WASH  DC  20314 

HQ  AF1SC/SEV  1 

Norton  AFB,  CA  92409 

OUTSIDE  AGENCIES 

Battelle,  Columbus  Laboratories 

Attn:  TACTEC  1 

505  King  Avenue 
Columbus,  OH  43201 

Director  of  Toxicology  1 

National  Research  Council 
2101  Constitution  Ave,  NW 
Washington,  DC  20418 

ADDITIONAL  ADDRESSEES 


Commander/Director 
Chemical  Systems  Laboratory 


Attn: 

DRDAR-CLG 

1 

Attn: 

DRDAR-CLL 

1 

Attn: 

DRDAR-CLC 

1 

Attn: 

DRDAR-CLN-S 

1 

Attn: 

DRDAR-CLW-P 

1 

Attn: 

D RD AR-CLB-PS(Mr.  Vervier) 

1 

Attn: 

Aberdeen 

DRDAR-CLB-PS(Dr.  Stuebing) 

Proving  Ground,  MD  21010 

1 

Names  Copies 

ADDITIONAL  ADDRESSEES 

Commander 

USEUCOM 

Attn:  ECJ-O/LTC  James  H.  Alley  1 

APO,  NY  09128 

US  Public  Health  Service  1 

Room  17A-46  (CPT  OsherofO 
5600  Fishers  Lane 
Rockville,  MD  20856 


Commander 

US  Army  Environment  Hygiene  Agency 

Attn:  Librarian,  Bldg  2100  1 

Aberdeen  Proving  Ground,  MD  21010 

Commander 
DARCOM,  ST1TEUR 

Attn:  DRXST-ST1  1 

Box  48,  APO  New  York  09710 

Commander 

US  Army  Science  & Technology  Center-Far  East  Office  1 
APO  San  Francisco  96328 

HQDA  DASG-RDZ  (SGRD-PL)  1 

WASH  DC  20314 

Project  Manager  for  Smoke/Obscurants 

Attn:  DRCPM-SMK  2 

Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Materiel  Systems  Analysis  Activity 

Attn:  DRXSY-D  (Dr.  Fallin)  1 

Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Ballistic  Research  Laboratories 

Attn.  DRXBR-DL(Mr.  T.  Finnerty)  1 

Attn:  DRXBR-P(Mr.  N.  Gerri)  1 

Attn:  DRDAR-BLB(Mr.  R.  Reitz)  1 

Attn:  DRDAR-BLB(Mr.  A.  LaGrange)  1 

Aberdeen  Proving  Ground,  MD  21005 

CDR,  APG 

Attn:  STEAP-TL  1 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Ordnance  and  Chemical  Center  and  School 

Attn:  ATSL-CD-MS  1 

Attn:  ATSL-CD-MS(Dr.  T.  Welsh)  1 

Aberdeen  Proving  Ground,  MD  21005 

Office  of  the  Director 

Defense  Research  and  Engineering 

Attn:  Dr.  T.C.  Walsh,  Rm  3D-1079  1 

Washington,  DC  20310 
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DISTRIBUTION  LIST  6 (Contd) 


Names  Copies 

Institute  for  Defense  Analysis 
400  Army-Navy  Drive 

Attn:  L.  Biberman  I 

Attn:  R.E.  Roberts  1 

Arlington,  VA  22202 

Commander 

US  Army  Research  Office  - Durham  I 

Box  CM,  Duke  Station 
Durham,  NC  27706 

HQ  (DAMA-ARZ,  Dr.  Verderame)  1 

HQDA  (DAMA-CSM-CM)  1 

WASH,  DC  20310 

Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Attn:  DRCDE-DM  1 

Attn:  DRCMT  1 

Attn:  DRCSF-S  1 

Attn:  DRCDL/Mr.  N.  Klein  1 

Attn:  DRCBI/COL  Gearin  1 

Attn:  DRCDMD-ST/Mr.  T.  Shirata)  1 

5001  Eisenhower  Ave. 

Alexandria,  VA  22333 


Commander 

US  Army  Foreign  Science  & Technology  Center 

Attn:  DRXST-ISI  1 

Attn:  DRXST-CE/Mr.  V.  Rague  1 

220  Seventh  St.,  NE 
Charlottesville,  VA  22901 


Commander 

US  Army  Missile  Command 

Attn:  Chief,  Documents  1 

Attn:  DRDMI-CGA/Dr.  B.  Fowler  1 

Attn:  DRDMI-TE/Mr.  H.  Anderson  1 

Attn:  DRDMI-KL/Dr.  W.  Wharton  1 

Redstone  Arsenal,  AL  35809 

Advanced  Research  Projects  Agency  1 

1400  Wilson  Boulevard 
Arlington,  VA  22209 


Commander 

Harry  Diamond  Laboratories 

Attn:  DRXDO-RDC(Mr.  D.  Giglio)  1 

2800  Powder  Mill  Road 
Adelphi,  MD  20783 

Chief,  Office  of  Missile  Electronic  Warfare 
US  Army  Electronic  Warfare  Laboratory 

Attn:  DRSEL-WLM-SE(Mr.  K.  Larson)  1 

White  Sands  Missile  Range,  NM  88002 

US  Army  Mobility  Equipment  Research  and 
Development  Center 

Attn:  Code/DROME-RT(Mr.  O.F.  Kezer)  1 

Fort  Belvoir,  VA  22060 


Names  Copies 

Commander 

Atmospheric  Sciences  Laboratory 

Attn:  DRSEL-BR-AS-P  1 

Attn:  DRSEL-BR-MS-A(Dr.  R.  Gomez)  1 

Attn:  DRSEL-BL-AS-DP(Mr.  J.  Lindberg)  1 

Attn:  DRSEL-BL-SY(Mr.  F.  Horning)  1 

White  Sands  Missile  Range,  NM  88002 

Director 

Night  Vision  Laboratories 

Attn:  DRSEL-NV-VI(Mr.  R.  Moulton)  I 

Attn:  DRSEL-NV-VI(Mr.  R.  Bergemann)  1 

Fort  Belvoir,  VA  23651 


Commander 

US  Army  Electronics  Command 

Attn:  DRSEL-CT-LG(Dr.  R.G.  Rohde)  1 

Attn:  DRSEL-CT-l(Dr.  R.G.  Buser)  1 

Attn:  DRSEL-WL-S(Mr.  J.  Charlton)  1 

Fort  Monmouth,  NJ  07703 

Commander 

Dugway  Proving  Ground 

Attn:  STEDP-PO  1 

Attn:  Technical  Library,  Docu  Sec  1 

Attn:  STEDP-MT-DA-E  1 

Attn:  STEDP-MT  (Dr.  L.  Salamon)  1 

Dugway,  UT  84022 


Commander 

US  Army  Logistics  Center 

Attn:  ATCL-MM  1 

Fort  Lee,  VA  23801 

Commander 
HQ,  USA  TRADOC 

Attn:  ATCD-TEC  (Dr.  M.  Pastel)  1 

Fort  MOnroe,  VA  23651 

Commander 

Naval  Surface  Weapons  Center 

Attn:  Tech  Lib  & Info  Svcs  Br  1 

White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Commander 

Naval  Intelligence  Support  Center  1 

4301  Suitland  Road 
Washington,  DC  20390 

Commander 

Naval  Surface  Weapons  center 
Dahlgren  Laboratory 

Attn:  DX-21  1 

Dahlgren,  VA  22448 

Commander 

Armament  Development  & Test  Center 

Attn:  DLOSL  (Technical  Library)  1 

Eglin  AFB,  FL  32542 
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DISTRIBUTION  LIST  6 (Contd) 


Names 

Copies 

Names 

Copies 

Commander 

Naval  Weapons  Center 

Attn:  Code  3311  (Dr.  R.  Bird) 

Attn:  Code  382  (Dr.  P.  St.  Amand) 

Attn:  Code  3822  (Dr.  Hindman) 

China  Lake,  CA  93555 

1 

1 

1 

Commander 

Naval  Research  Laboratory 

Attn:  Code  5709  (Mr.  W.E.  Howell) 

4555  Overlook  Avenue,  SW 

Washington,  DC  20375 

1 

Commanding  Officer 

Naval  Weapons  Support  Center 

Attn:  Code  5041  (Mr.  D.  Johnson) 

1 

US  Army  Materiel  Systems  Analysis  Activity 

Attn:  DRXSY-MP 

Aberdeen  Proving  Ground,  MD  21005 

1 

Crane,  IN  47522 
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